Phonon and crystal field excitations in geometrically frustrated rare earth titanates 
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The phonon and crystal field excitations in several rare earth titanate pyrochlores are investigated. 
Magnetic measurements on single crystals of Gd2Ti207, Tb2Ti207, Dy2Ti207 and Ho2Ti207 are 
used for characterization, while Raman spectroscopy and terahertz time domain spectroscopy are 
employed to probe the excitations of the materials. The lattice excitations are found to be analogous 
across the compounds over the whole temperature range investigated (295-4 K). The resulting 
full phononic characterization of the R2Ti207 pyrochlore structure is then used to identify crystal 
field excitations observed in the materials. Several crystal field excitations have been observed in 
Tb2Ti207 in Raman spectroscopy for the first time, among which all of the previously reported 
excitations. The presence of additional crystal field excitations, however, suggests the presence of 
two inequivalent Tb^^ sites in the low temperature structure. Furthermore, the crystal field level at 
approximately 13 cm~^ is found to be both Raman and dipole active, indicating broken inversion 
symmetry in the system and thus undermining its current symmetry interpretation. In addition, 
evidence is found for a significant crystal field-phonon coupling in Tb2Ti207. These findings call for 
a careful reassessment of the low temperature structure of Tb2Ti207, which may serve to improve 
its theoretical understanding. 



PACS numbers: 63.20.dd, 71.70.Ch, 75.50.Lk, 75.30.Cr 
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I. INTRODUCTION 

The term 'geometrical frustration 'i^i^ applies to a sys- 
tem, when it is unable to simultaneously minimize all 
of its magnetic exchange interactions, solely due to its 
geometry. Magnetically interacting spins residing on 
such lattices are unable to order into a unique magnetic 
ground state due to the competing magnetic interactions 
between different lattice sites. Instead of selecting a sin- 
gle, unique magnetic ground state at low temperatures, 
a pure magnetically frustrated system has a macroscop- 
ically degenerate ground state. In real systems however, 
any secondary, smaller term (arising from single-ion or 
exchange anisotropy, further neighbor interactions, dipo- 
lar interactions, small lattice distortions or a magnetic 
field, for example) in the system's Hamiltonian can fa- 
vor certain magnetic ground states at very low tempera- 
tures, thereby (partially) lifting this peculiar degeneracy. 
In this fact lies the origin of the vast richness and diver- 
sity of the low temperature magnetic behavior of different 
frustrated systems in naturei'^i^'^ii. 

Geometries suitable to exhibit frustration typically 
consist of infinite networks of triangles or tetrahedra, 
which share one or more lattice sites. One of the most 
common structures known to be able to induce mag- 
netic frustration is that of the pyrochlores, A2B2O7, 
where both the A'^^ ions (rare earth element, coordi- 
nated to 8 O atoms) and the B^+ ions (transition metal 
element, coordinated to 6 O atoms) reside on a lattice of 
corner-sharing tetrahedra, known as the pyrochlore lat- 
tice. Thus, if either A'^+ or B"'+ is a magnetic species. 



frustration may occur due to competing interactions. A 
subclass of the pyrochlores is formed by the rare earth 
titanate family, R2Ti207, where the R^+ ion is the only 
(para-)magnetic species, since Ti*^"*" is diamagnetic {3dP). 
For the pyrochlore lattice, both theor}*2ii2iiiS and Monte 
Carlo simulationsiSiiii predict a 'collective paramagnetic' 
ground state, or the lack of long range magnetic order- 
ing, for classical Heisenberg spins at finite temperature. 
The quantum Heisenberg spin (S = 1/2) model for the 
pyrochlore lattice also predicts a quantum disordered sys- 
tem at finite temperatures, a state often referred to as a 
'spin liquid'—. However, in reality the different pertur- 
bative terms in the corresponding Hamiltonian result in 
quite diverse low temperature magnetic behavior among 
the different rare earth titanates^, of which the Gd, Tb, 
Ho and Dy variants are studied here. 

The supposedly least complex case is that of gadolin- 
ium titanate, Gd2Ti207. The Gd"^^ ion has, in contrast 
to the Tb^+, Ho^+ and Dy^+ ions, a spin only ^87/2 
[L = 0) ground state, rendering the influence of crys- 
tal field levels and possible induced Ising-like anisotropy 
insignificant in Gd2Ti207. The experimentally deter- 
mined Gurie-Weiss temperature of Gd2Ti207 is ci —10 
^13^14,15 ^ indicating antiferromagnetic nearest neighbor 
interactions. Thus, Gd2Ti207 could be considered as 
an ideal realization of the frustrated Heisenberg anti- 
ferromagnet with dipolar interactions. Experimentally, 
Gd2Ti2 07 has been found to undergo a magnetic order- 
ing transition at ~ 1 K-i2.. However, this transition cor- 
responds to only partial ordering of the magnetic struc- 
ture, as only 3 spins per tetrahedron order—. In this 



partially ordered state the spins residing on the [111] 
planes of the crystal (which can be viewed as Kagome 
planes) are ordered in a 120° configuration, parallel to the 
Kagome plane, while the spins residing on the interstitial 
sites remain either statically or dynamically disorderedii. 
Subsequent experimental investigations revealed a sec- 
ond ordering transition at ~ 0.7 K, corresponding to the 
partial ordering of the interstitial disordered spinsi^iii 
which, however, do remain (partially) dynamic down to 
20mK .^^'^^ Despite of Gd2Ti207 supposedly being well 
approximated by the Heisenberg antiferromagnet with 
dipolar interactions, theoretical justification for this com- 
plex magnetic behavior remains difficul t^^i^^i^° . 

In Tb2Ti207, the dominant interactions are antifer- 
romagnetic, as indicated by the experimentally deter- 
mined Curie- Weiss temperature, Ocw — —19 K^i. A 
study of the diluted compound (Tbo.o2Yo.g8)2Ti207 re- 
vealed that the contribution to 6cw due to exchange 
and dipolar interactions is ~ —13 K, comparable to 
the 0cw-value found in Gd2Ti207^. Despite the en- 
ergy scale of these interactions, the Tb^+ moments do 
not show long range magnetic order down to as low as 
50 niK, making it the system closest to a real 3D spin 
liquid to date^iiS^. However, crystal field (CF) calcula- 
tions indicate a ground state doublet and Ising-like easy 
axis anisotropy for the {J Fq) Tb^"*" magnetic moments 
along their local <111> directions (the direction towards 
the center of the tetrahedron the particular atom is in), 
which would dramatically reduce the degree of frustra- 
tion in the system^S'^i^'^^i^. Theoretical models taking 
this anisotropy into account predict magnetic ordering 
temperatures of about 1 K^^'^^ . Subsequent theoreti- 
cal work suggests that the magnetic moment anisotropy 
is more isotropic than Ising-like, which could suppress 
magnetic ordering22i. Recently, Tb2Ti207 was argued 
to be in a quantum mechanically fluctuating 'spin ice' 
state^i^. Virtual quantum mechanical CF excitations 
(the first excited CF doublet is separated by only ~ 13 
cm~^ from the ground state double t^^'^^i^^ ) are proposed 
to rescale the effective theoretical model from the un- 
frustrated Ising antiferromagnet to a frustrated resonat- 
ing spin ice model. Nevertheless, the experimentally ob- 
served lack of magnetic ordering down to the millikelvin 
range and the true magnetic ground state in Tb2Ti207 
remain as of yet enigmatio^i^^i^i^i. 

Illustrating the diversity in magnetic behavior due to 
the subtle differences in the rare earth species of the 
titanates, the situation in both Dy2Ti20 ,^^i'^^'^^ and 



Ho2Ti207 
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is again different. The R'^"'" ions 



in these compounds have a ®i?i5/2 (Dy'^+) and a ^/g 
(Ho'^"'') ground state, respectively, with corresponding 
free ion magnetic moments ^ = 10.65 ^b (Dy'^"'') and 
/I — 10.61 ^b{)^o^^)- These systems were first thought 
to have weak ferromagnetic nearest neighbor exchange 
interactions, as indicated by the small positive values 
of Ocw, - 2 K and ~ 1 K for Dy2Ti207 and 
Ho2Ti207, respectively^. More recently, however, the 
nearest neighbor exchange interactions in Dy2Ti207 and 



Ho2Ti207 were argued to be aniiferromagnetio^. The 
effective ferromagnetic interaction between the spins in 
fact is shown to be due to the dominant ferromagnetic 
long-range magnetic dipole-dipole interactions^ -'^^ . The 
R'^+-ions in both Dy2Ti207 and Ho2Ti207 are well de- 
scribed by a well separated Ising doublet (first excited 
states are at ~ 266 and ^ 165 cm~^ ,respectively2i) 
with a strong single ion anisotropy along the local < 
111 > directions. Unlike for antiferromagnetically in- 
teracting spins with local < 111 > Ising anisotropy, 
/erromagnetically interacting Ising spins on a pyrochlore 
lattice should be highly frustrated^i^. As Anderson 
already pointed out half a century agc4^, the result- 
ing model is analogous to Pauling's ice model^, which 
earned both Dy2Ti207 and IIo2Ti207 the title of 'spin 
ice' compound^i^i^i^. Although numerical simulations 
predict long range order at low temperatures for this 
model^, experimental studies report no transition to a 
long range ordered state for either Dy2Ti20 |^^i'^^'^^i^° or 
Ho2Ti207^i^iS, down to as low as 50 mK. 

As is apparent from above considerations, the low tem- 
perature magnetic behavior of the rare earth titanates is 
dictated by the smallest of details in the structure and 
interactions of the material. Therefore, a comprehensive 
experimental study of the structural, crystal field and 
magnetic properties of these systems may serve to clar- 
ify unanswered questions in their understanding. In this 
paper, dc magnetic susceptibility measurements, polar- 
ized Raman scattering experiments and terahertz time 
domain spectroscopy on aforementioned members of the 
rare earth titanates family are employed to gain more in- 
sight into the details that drive them towards such diverse 
behavior. Raman scattering allows for simultaneous in- 
vestigation of structural and crystal field (CF) properties 
through the observation of both phononic and CF exci- 
tations, while the comparison between the various mem- 
bers helps identify the nature of the different excitations 
observed. 



II. EXPERIMENTAL 
A. Sample Preparation 

Polycrystalline samples of R2Ti207 (where R = Gd, 
Tb, Dy, Ho) were synthesized by firing stoichiometric 
amounts of high purity (> 99.9%) Ti02 and the appropri- 
ate rare earth oxide (Gd203, Tb407, Dy203 or H02O3, 
respectively), in air, for several days with intermittent 
grindings. The resulting polycrystalline powder was sub- 
sequently prepared for single crystal growth, using the 
method described by Gardner, Gaulin and Paul^. The 
following single crystal growth (also done as described in 
ref. 49) using the floating zone technique yielded large, 
high quality single crystals of all of the R2Ti207 variants. 
Discs (~ 1 mm thickness) with a,b-plane surfaces were 
cut from oriented single crystals and subsequently pol- 
ished, in order to optimize scattering experiments. The 



Tb2Ti207 sample used in Raman experiments was subse- 
quently polished down to ~ 250 /ini thickness to facilitate 
THz transmission measurements. 



B. Instrumentation 

X-ray Laue diffraction, using a Philips PW 17f0 
diffractometer equipped with a Polariod XR-7 system, 
was employed to orient the single crystal samples of 
Gd2Ti207 and Tb2Ti207 for the polarized Raman spec- 
troscopy experiments, while simultaneously confirming 
the single crystallinity of the samples. The Dy2Ti2 07 
and Ho2Ti207 single crystals were oriented using an En- 
raf Nonius CAD4 diffractometer. 

The magnetic susceptibilities of the obtained rare 
earth titanates were measured using the Quantum De- 
sign MPMS-5 SQUID magnetometer of the 'Laboratoire 
de Physico-Chimie de I'Etat Sohde' (LPCES), CNRS, 
UMR8182 at the 'Universitc Paris-Sud in Orsay, France. 
The R2Ti207 samples, about 100 mg of single crystal 
(in the form of discs of approximately 4 mm diameter 
and 1 mm thickness), were placed in cylindrical plastic 
tubes and locked in position. Next, the samples were 
zero-field-cooled down to 1.7 K, after which the magne- 
tization of the sample was measured as a function of the 
temperature in an applied magnetic field of 100 Oe, while 
warming the sample. 

Polarization controlled, inelastic light scattering exper- 
iments were performed on all oriented R2Ti207 samples. 
The experiments were performed in a 180° backscattering 
configuration, using a triple grating micro-Raman spec- 
trometer (T64000 Jobin-Yvon), consisting of a double 
grating monochromator (acting as a spectral filter) and 
a polychromator which disperses the scattered light onto 
a liquid nitrogen cooled CCD detector. The frequency 
resolution was better than 2 cm~^ for the frequency re- 
gion considered. The samples were placed in a liquid he- 
lium cooled optical flow-cryostat (Oxford Instruments). 
The temperature was stabilized with an accuracy of 0.1 
K in the whole range of measured temperatures (from 
2.5 to 295 K). The 532.6 nm (frequency doubled) out- 
put of a Nd:YV04 laser was focused on the Gd2Ti207, 
Tb2Ti207 and Dy2Ti207 samples using a 50x microscope 
objective and used as excitation source in the scattering 
experiments. A Krypton laser (676.4 nm) was used as 
the excitation source for the scattering experiments on 
the Ho2Ti207 sample, since 532.6 nm excitation (reso- 
nant at low temperatures in case of IIo2Ti207) results 
in fluorescence dominating the inelastic scattering spec- 
trum in the 5-800 cni^^ spectral range. The power den- 
sity on the samples was of the order of 50 /j,W/^m^ in 
all cases. The polarization was controlled both on the 
incoming and outgoing beam. Parallel (||) and perpen- 
dicular (_L) measurements on Gd2Ti207 and Tb2Ti207 
were performed along crystallographic axes of the a,b sur- 
face of the samples, Porto notations c(aa)c and c(ab)c, 
respectively. Unfortunately, the orientation of the a and 



b axes in the (a,b) plane of the Dy2Ti207 and IIo2Ti207 
surfaces with respect to the light polarizations was not 
known. Analogous Porto notations are c(xx)c (||) and 
c(xy)c (_L), respectively, where x is a direction in the 
(a,b) plane of the sample making an undetermined an- 
gle a with the a axis, while y, in the same a,b plane of 
the crystal, is perpendicular to the x direction. Raman 
spectra were fitted with lorentzian lineshapes to extract 
mode parameters. 

Terahertz time domain spectroscopy (THz TDS)^ was 
performed on Tb2Ti207 using a home made setup simi- 
lar to those described elsewhere^^*^. THz pulses (pulse 
duration of several ps, frequency range 0.3-2.5 THz) were 
generated through a difference frequency generation pro- 
cess in a ZnTc single crystal upon pulsed excitation (120 
fs. 800 nm) by an amplified Ti:sapphire system. The 
magnitude of the time dependent electric field transmit- 
ted through the sample (w.r.t. that transmitted through 
vacuum) was measured at various temperatures, through 
electro-optic sampling in a second ZnTe single crystal, us- 
ing 800 nm pulses of approximately 120 fs. The sample 
used was a thin slice of single crystalline Tb2Ti207 (the 
same sample as used in the Raman experiments), which 
was mounted on a copper plate with an aperture (0 2 
mm) and placed in a liquid helium cooled optical flow- 
cryostat (Oxford Instruments). The polarization of the 
THz radiation was parallel to the crystallographic a-axis. 



III. RESULTS AND DISCUSSION 

A. Magnetic Measurements 

The magnetic susceptibility x, defined as the ratio of 
the magnetization of the sample to the applied magnetic 
field, of all the rare earth titanate samples was measured 
in a 100 Oe applied magnetic field. Since the samples 
used were plate-like discs, the data have been corrected 
by a demagnetization factor as calculated for flat, cylin- 
drical plates^. Fig. [T] shows the inverse molar suscepti- 
bilities of all samples in the low temperature regime. For 
each sample, the data was fitted to a Curie- Weiss form for 
the molar magnetic susceptibility of an antiferromagnet: 
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where C is the Curie constant in CGS units (C ~ 
NAfJ-^H^/Skij ~ ^^/8, in [emu-K-mol~^]), 9 is the ex- 
pected transition temperature (giving an indication of 
the sign of the magnetic interactions) and i? is a tem- 
perature independent Van Vleck contribution to the sus- 
ceptibility. The model was fitted to high temperature 
experimental data (100 K and up, where the demagne- 
tization correction is of the order of 1 %) and linear re- 
gression analysis of Xm i yielded the experimental values 
for 9i, fii and B^. These are tabulated below in table 1, 
together with several values reported in literature. 
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FIG. 1: Inverse molar magnetic susceptibilities (open circles) of all R2Ti207 samples and their corresponding Curie- Weiss fits 
(solid lines, eq. [T| versus temperature. The Curie- Weiss law, which is fitted to the higher temperature data (> 100 K) is 
obeyed down to low temperatures. 



TABLE I: Experimentally determined and literature values 
for 6i, fj,i and Bi for the different R2Ti2 07 compounds. The 
numbers in between brackets give the accuracy in the last 
decimal of the values. 



Compound i (source) 


e.{K) 


tJ-ii/J-B) B, 


{emu ■ mol 


") 


Gd2Ti2 07 (exp.) 


-11(2) 


7.95(9) 


0" 




Gd2Ti207 (lit.y^ 


-9,6 


7.7 


- 




Gd2Ti207 (lit.)^ 


-11,7 


7.8 


- 




Gd2Ti207 (ht.)^ 


-8,95(6) 


7.224(3) 


6.0 • 10 


-4 


Tb2Ti207 (exp.) 


-18(2) 


9.67(7) 


3.4(9) ■ 


10- 


Tb2Ti207 (te.)2i-2a 


-19 


9.6 


- 




Dy2Ti2 07 (exp.) 


2(2) 


6 


b 




Dy2Ti207 (te.)^ 


-1,0 


9.590(6) 


1.1(1). 


10" 


Ho2Ti207 (exp.) 


6(2) 


_h 


_b 




H02Ti207 (Itt.)^ 


-2,0 


9.15(3) 


1.2(1). 


10- 



"When fitting Xm with a nonzero B term, it becomes negligibly 
small. The best fit is obtained with B=0. 

''Since dipolar interactions are dominant in this case, no reliable 
values were extracted (see text). 



In general, the experimentally obtained data compare 
(where possible) favorable to the various values reported 
in literature (table |T| . The extracted parameters for 
Dy2Ti207 and Ho2Ti207 do slightly deviate from litera- 
ture values, presumably due to the estimation of the de- 
magnetization factor (fits to the present uncorrected data 
yield 9 values of approximately 1 and 2 K, respectively). 
The experimentally determined paramagnetic moments 
obtained for Gd2Ti207 and Tb2Ti207 are also in excel- 
lent agreement with the corresponding free ion values, 
which are // = 7.94 /i^ and /i = 9.72 jib for the Gd'^+ 



(^5*7/2) and Tb'^'*" ij Fq) free ions, respectively. The large 
negative Curie-Weiss temperatures for Gd2Ti207 and 
Tb2Ti207 indicate antiferromagnetic exchange coupling. 
In contrast, the small, positive 9 values for Ho2Ti207 and 
Dy2Ti2 07 initially led to the assumption of weak ferro- 
magnetic exchange interactions between nearest neigh- 
bor Dy'^"'' and IIo"^+ ions^. As stated above, however, 
since the Ho'^+ and Dy'^+ ions have a large magnetic mo- 
ment (free ion values are /i = 10.607 ^b (^-^s ground 
state) and ^ = 10.646 /zb (^^15/2 ground state), re- 
spectively), the dipolar interactions between neighboring 
R'^+ ions are dominating the effective nearest neighbor 
(n.n.) interactions. The n.n. exchange interactions are 
in fact antzferromagnetio^, while the dominant dipolar 
n.n. interactions are of ferromagnetic nature^i^. Con- 
sequently, the effective n.n. interactions are slightly fer- 
romagnetic, resulting in the positive 9 values. Another 
consequence of the dipolar interactions and being domi- 
nant is the fact that extracting the real values of /i and 
B from the inverse susceptibility curves becomes non- 
trivial, since more elaborate models taking the dipolar 
interaction into account are needed. 



B. Raman Spectroscopy 

1. Room temperature spectra 

To determine the Raman active vibrations of the sin- 
gle crystals, group theory analysis was employed. This 
predicts that, for the cubic rare earth titanate structure 
(R2Ti2 07) of space group FdSm (OfJ), the sublattices 



of the unit cell span the following irreducible representa- 
tions: 



16(c)-site: Ti^+ — sublattice = 

16((i)-site: R^+ — sublattice = 

48(/)-site: 0(1) - sublattice = 

8(a)-site: 0(2) — sublattice = 
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This makes the following decomposition into zone cen- 
ter normal modes (excluding the Fi„ acoustic mode): 



T = Ai„ + E„ + 2Fu, + 4F2„ + 3A2„ + ZEu + 7Fi„ + AF, 



2m 



Of these normal modes, only the Aig, Eg and the 4 F2g 
modes are Raman active. The 7 Fiu modes are infrared 
active and the remaining modes are optically inactive. 
The symmetry coordinates for the optically active normal 
modes are given by Gupta et al^. Based on the symme- 
tries of the Raman active modes, the Aig and Eg modes 
are expected to be observed in parallel polarization (||) 
spectra, while the F2g modes are expected in the perpen- 
dicular polarization (_L) spectra. Surprisingly, the room 
temperature spectra of R2Ti207 show all the Raman ac- 
tive modes in both polarizations. Angle dependent Ra- 
man measurements (varying the angle between the in- 
coming light polarization and the b-axis from —45° to 
45°) reveal no orientation in which the theoretical selec- 
tion rules are fully obeyed. Polycrystallinity and/or crys- 
talline disorder is however, not believed to be the cause, 
since other techniques (X-ray Laue diffraction, magnetic 
measurements) reveal no such signs. Also, to our best 
knowledge, no Raman spectrum of any of these rare earth 
titanates fulfilling their phononic selection rules has been 
published. Although no full compliance of phononic se- 
lection rules was observed, the nature of the Raman ac- 
tive modes was clearly identified, through their angular 
dependence. Their assignment is indicated in fig. [5] as 
well as in table |TT1 The figure depicts the room tempera- 
ture (RT), parallel polarization spectra of the rare earth 
titanates and the corresponding modes. The assignment 
of the mode symmetries is based on comparison with as- 



signments in previous works on thes o^^'^^i^^ i 
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and 
relate d^^i^^'^'^1^^'^^1^^ compounds (see table |IT]), on afore- 
mentioned measurements (not shown) of the angular de- 
pendence of mode intensities and on the temperature de- 
pendence of the modes (vide infra). 

Comparison between the RT Raman spectra yields the 
conclusion that the nature of the R^+ ion has only slight 
infiuencc on the Raman active vibrational modes. This 
is not surprising, since in all the Raman active modes 
only the oxygen atoms are displaced^i^. Consequently, 
there is no obvious systematic variation of the phonon fre- 
quencies with the mass of the respective rare earth ions, 
as has been noted before for the rare earth titanates^, 
hafnatea^, manganatea^ and stannatea^. The assign- 
ment of the R2Ti2 07 modes in literature has been mostly 
consistent (see Table HI)) . yet there are a few debated de- 
tails. There is general agreement on the nature of the 



modes at ~ 210 cm ^ (F2g, 0(2)-sublattice mode^), ~ 
519 cm-i {Aig, R-0 stretching mode^i^) and ~ 556 
cm^^ (-^"29, 0(l)-sublattice mode^). Temperature and 
angle dependent Raman measurements show the band 
around ~ 315 cni^^ to consist of two modes, a F2g mode 
around 310 cm^^ (0-R-O bending mod a^^'^° ) and an Eg 
mode around 327 cm~^ (O(l)-sublattice mode^), as rec- 
ognized by Saha et al^ and Vandenborre et al^. Earlier 
works either interchanged the mode assignment within 
this band^i^ or ascribed the whole band to only one of 
these modes^i^. However, our temperature and angle 
dependent measurements confirm the assignment made 
by Saha et al. The last expected phonon, an F2g mode, 
has been either not accounted for—, combined with the 
Aig mode in one band^*^ or ascribed to low intensity 
peaks around 105^, 450^ or 680^ cm~^ in previous 
works. Here, it is ascribed to a broad, low intensity 
mode around ~ 260 cm"-'^. This mode is not clearly re- 
solved in the RT spectra due to the fact that it overlaps 
largely with the neighboring, strong i^2g modes at ~ 210 
cm~^ and ~ 309 cm~^. Fitting with those two peaks 
only, however, does not adequately reproduce the exper- 
imental spectral shape in the the 200-300 cm~^ window. 
Additionally, as the temperature is lowered, the phonon 
modes sharpen and the existence of this excitation be- 
comes obvious in the spectra. Worth noting are also the 
two anomalous modes in Tb2Ti207 at ~ 303 (^29) and ~ 
313 cm^^ {Eg), which have lower frequencies and wider 
lineshapes compared to their counterparts in the other, 
isostructural rare earth titanates. 

Next to the expected Raman active vibrations, the 
spectra in this work show some very weak scattering in- 
tensity at low wavenumbers (first two rows in table [TT)l . 
which has been reported beforo^i^. Vandenborre et al^ 
were unable to account for this intensity in their calcula- 
tions, while Mori et al^ offer the plausible assignment to 
trace R2O3 in the system. The latter assignment is also 
tentatively adopted here. Mori et al. also suggested the 
'missing' F2g mode may be responsible for some of this 
low wavenumber intensity. Additionally, a weak mode is 
observed at 450 cm~^ in most R2Ti207 compounds, as 
was also seen before. Zhang et al^ and Hess et al^ 
ascribed the 'missing' F2g mode to this feature. Alterna- 
tively, Mori et al^ interpreted it as being due to trace 
amounts of starting compound Ti02, which is known^ to 
have a phonon at 447 cm~^. The true origin of this mode 
is at present unclear. Finally, there is some low intensity 
scattering intensity at higher wavenumbers, around 700 
cm~^ (last two rows of table HI]) . This intensity has been 
observed befor a^^i^^i^°i^^i^^ and is ascribed to forbidden 
IR modes made active by slight, local non-stoichiometry 
in the systemi^i^ 



2. Temperature dependence 

Raman spectra of the R2Ti207 crystals were recorded 
at temperatures ranging from RT to 4 K. Figure[3]depicts 
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FIG. 2: (Color online) Room temperature Raman spectra of the R2Ti207 crystals in parallel polarization configuration, 
assignment of the Raman active modes to the observed peaks is indicated by the colored lorentzian lineshapes. 
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the evolution of both the parallel (||) and perpendicular 
(_L) polarization Raman spectra of the R2Ti207 crystals 
vifith decreasing temperature. 

Going down in temperature several spectral changes 
occur in the Raman spectra. The evolution of the 
phononic excitations with temperature are very similar in 
the different R2Ti207 lattices. Again, this is not surpris- 
ing, since only oxygen atoms are displaced in the Raman 
active phonons^i^. Firstly, the lowest frequency phonon 
{F2g, Lo ~ 210 cm"-'^ at RT) shows a strong softening 
(lo ^ 170 cm~^ at 4 K) and sharpening with decreas- 
ing temperature, revealing the previously unresolved F2g 
mode {ojRT — 260 cm~^), which also softens {uj4k — 190 
cm^^), but doesn't show a strong narrowing. Secondly, 
the sharpening of both the ~ 309 cm^^ F2g and the ~ 324 
cm""'^ Eg mode decreases their spectral overlap, clearly 
justifying the two-mode interpretation of the ~ 315 cm~^ 
band at RT. Additionally, both modes show a slight soft- 
ening upon cooldown. Also the Aig phonon {ujrt — 519 
cm""'^) shows the familiar softening {u!4k — 511 cm^^) 
and sharpening trend on cooling. Finally, due to its large 
width and low intensity, describing the temperature evo- 
lution of the highest frequency i^2g phonon proves rather 
difficult, though it seems to soften slightly. Compari- 
son of the R2Ti207 spectra in fig. [3] yields the observa- 
tion that the anomalous phonons (at ~ 303 {F2g) and 
~ 313 cm^^ {Eg)) ill Tb2Ti207 remain wide throughout 
the temperature range, in contrast to the corresponding 
modes in the other titanates. Additionally, these modes 
are shifting in opposite directions in Tb2Ti207 only: the 
F2g mode softens {uj4k — 295 cm~^) while the Eg mode 
considerably hardens {uj4k — 335 cm~^). An explana- 
tion for this anomalous behavior could be coupling of 



these phonons to low frequency crystal field excitations 
of the Th'^^-ions (vide infra). 



3. Crystal field excitations 

Aside from the phonons in the Raman spectra of 
R2Ti207, several spectra also show crystal field (CF) ex- 
citations of the R'^"'" ions at low temperatures, as shown 
in fig. m The CF level splitting in the different rare earth 
ions depends on their electronic configuration and their 
local surroundings. In the R2Ti207 family, the simplest 
case is that of the Gd'^^ ion (4/^), which has a spin only 
^5*7/2 {L = 0) ground state, resulting in the absence of a 
level splitting due to the local crystal field. Consequently, 
the Raman spectrum of Gd2Ti207 shows no CF exci- 
tations, making it a suitable 'template' of the R2Ti207 
Raman spectrum with lattice excitations only. Combined 
with the strong correspondence of the phononic excita- 
tions in the R2Ti207 spectra, it allows for quick identifi- 
cation of CF modes in the other compounds. 

More complicated is the CF level splitting in the 
Tb'^^ ion (4/^), which has a ^Fg ground state. Sev- 
eral studies calculating the crystal field for the Tb^^^-ion 
in Tb2Ti20 , 22i25, 26^27 ^ based on inelastic neutron scat- 
tering rcsult a^-'^'^^i^^'^'^i^^''^" , yielded slightly differing en- 
ergy level schemes for the lowest crystal field levels of 
Tb'^+ in Tb2Ti207 as is schematically depicted in fig. 
[5l Here, solid lines depict CF levels observed experi- 
mentally, while dotted lines indicate CF levels obtained 
through CF calculations. Shown in fig. [SI which is a 
zoom-in on the low-wavcnumbcr region of the 4 K per- 
pendicular polarization spectrum of Tb2Ti207, are the 
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FIG. 3; Temperature dependence of the R2Ti207 parallel (||) and perpendicular {!.) Raman spectra. The spectra show several 
changes when going down in temperature (see text), most of which (the phononic excitations) are analogous for all R2Ti2 07 
lattices. In addition, at lower temperatures various crystal field excitations appear in the spectra of some of the rare earth 
titanates (See fig. [J] Spectra are normalized to the integrated intensity of the depicted frequency window. 



CF excitations that are observed using inelastic light 
scattering. As is also clear from fig. [5l all CF excita- 
tions previously observed using inelastic neutron scat- 
tering are also observed here. Furthermore, additional 
lo'w-lying excitations can be seen, "which are easily iden- 
tified as CF levels, by comparison -with the 'lattice-only 



template' spectrum of Gd2Ti207. 

The excitations from the crystal field ground state to 
the higher crystal field levels, of approximate calculated 
values 13, 60, 83 and 118 cm~^ (sec fig. [S]), are all ob- 
served at very similar frequencies, at 12.9, 60.8, 83.1 and 
119.2 cm"-'^, respectively. The 60.8 cm~^ level has not 
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FIG. 4: (Color online) Low temperature spectra for all R2Ti207 as measured in perpendicular polarization configuration. The 
spectra are normalized on the total integrated intensity of the spectra. The insets show respective zoom-ins on the lower 
wavenumber regions of the spectra, where the excitations identified as CF levels of the rare earth ions are indicated by dotted 
blue lorentzian peaks. Gd2Ti207 shows only phononic modes and can be regarded as the lattice template. Tb2Ti207 and 
Dy2Ti207 show additional CF modes, while Ho2Ti207 shows again only lattice excitations. A fitted lorentzian centered at 
zero (Rayleigh line) is subtracted from the Tb2Ti207 data for clarity. 



been observed experimentally before, though it did fol- 
low from the CF calculation made by Gingras et ali^. 
Conversely, the ~ 135.2 cm~^ mode has been observed 
through inelastic neutron scattering, yet has not been 
accounted for in CF calculations. Although Gardner et 
al^ interpreted it as an optical phonon, this excitation 
is clearly identified as a CF excitation here, through the 
isostructural comparison with the other titanates. The 
latter assignment is also made by Mirebeau et al^. Ad- 
ditionally, new CF excitations are observed at 7.7 and 
47.2 cm~^. While the former is recognized as a CF ex- 
citation from the CF ground state to a new CF energy 
level, the 47.2 cm~^ could also be interpreted as an ex- 
citation from the excited CF level at 12.9 cm^^ to the 
higher excited CF level at 60.8 cm~^. While such an ex- 
citation is Raman active (see below), its occurrence at 4 
K seems unlikely, since at this temperature the excited 
CF level (12.9 cm~^ ~ 19 K) is not expected to be popu- 
lated enough to give rise to a measurable Raman signal. 
Additionally, were it an excitation from an excited state, 
its intensity would decrease upon cooling. Instead, the 
intensity of this mode steadily increases with decreasing 
temperature. 



Using Eg and Eg, A2g and A 



representation; 



,22.26 



Us 



irreducible 



for the ground state and excited CF 



levels, respectively, it can be confirmed that indeed all 
these ground state to excited state transitions are ex- 



pected to be Raman active. Through analogous consid- 
erations, the lowest excitation {Eg -^ Eg, 13 cm~^) is 
found to be symmetry forbidden in a direct dipole transi- 
tion. In this respect, terahertz time domain spectroscopy 
(THz TDS) in the range of 10 to 50 cm^^ (0.3 to 1.5 
THz) was employed to probe the absorption of Tb2Ti207 
at various temperatures. Using this technique it is pos- 
sible to extract complex optical quantities in the THz 
range through the direct measurement of the time trace 
of the transmitted THz pulse. The obtained curves for 
the real [eiiv)) and imaginary (£2(1^)) parts of the dielec- 
tric constant at various temperatures are plotted in fig. 
[T] The plots clearly show an absorption around 0.45 THz 
(15 cm"-"^) at low temperatures, corresponding to the 13 
cm"-*^ CF level also observed in the Raman spectra of 
Tb2Ti207. This is corroborated by the identical tem- 
perature dependence of the mode in both methods: its 
spectral signature decreases with increasing temperature, 
vanishing around 90 K. The fact that the 13 cm^^ CF 
level is observed in a direct dipole transition indicates 
that this level is not a true Eg -^ Eg transition, since 
such a dipole transition would be symmetry forbidden. 

Overall, the existence of the additional CF levels in 
Tb2Ti207 unaccounted for by CF calculations, suggests 
the presence of a second Tb^+ site in the structure, with 
an energy level scheme different from those reported pre- 
viously. Three additional CF levels are observed exper- 
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FIG. 5; Various CF energy level schemes of the Tb"^^ ion in Tb2Ti207. Solid lines depict CF levels that have been observed 
experimentally, while dotted lines represent CF levels as calculated using CF calculations. Numbers indicate CF level energies 
in cm~^ (THz) and symbols in italics indicate symmetries of corresponding levels. Solid arrows indicate all CF excitations 
observed in the low temperature Raman spectrum of Tb2Ti207 (see fig. |6]). 



imcntally, while four would naively be expected for a 
slightly differing Tb site in this low wavenumber region. 
A fourth new CF level might be unresolved due to the 
strong CF level at 83 cm""'^ or might simply be symmetry 
forbidden in a Raman transition. Moreover, the fact that 
the ~ 13 cm^^ CF level is simultaneously Raman and 
dipole active indicates the breaking of inversion symme- 
try in the system, which questions the validity of its cur- 
rent symmetry interpretation. Recently, the exact sym- 
metry of the Tb2Ti2 07 lattice has been extensively stud- 
ied. Han et ali^ performed neutron powder diffraction 
and x-ray absorption fine-structure experiments down to 
4.5 K, revealing a perfect pyrochlore lattice, within ex- 
perimental error. Ofer et ali^ find no static lattice dis- 
tortions on the timescalc of 0.1 /is, down to 70 mK. Most 
recently, however. Ruff et al^ found finite structural cor- 
relations at temperatures below 20 K, indicative of fluc- 
tuations above a very low temperature structural tran- 
sition. The present experimental results suggest these 
fluctuations may even induce a minute static disorder, 
resulting in the observed CF level diagram. 

In DyaTiaO? the Dy3+ ions (4/^) have a ^7Ii5/2 ground 



state. Crystal field calculations have been performed by 
Jana et al^, who deduced an energy level scheme con- 
sisting of eight Kramers' doublets, with a separation of 
~ 100 cm"-'^ of the first excited state. Malkin et ali^ and 
Rosenkranz et alJ^ estimated the first excited state gap 
to be > 200 cm^^ and ~ 266 cm^^, respectively. The 
Raman spectrum of Dy2Ti207 shows only one extra ex- 
citation compared to the 'lattice template' spectrum of 
Gd2Ti207, at an energy of ^ 287 cm~^. This excitation 
is tentatively ascribed to the first excited CF level, com- 
paring most favorably to the estimation of Rosenkranz et 
al.. 

For HoaTiaOr, ground state ^h (Ho3+,4/i0), several 
CF energy level schemes have been calculate d^^'^^i^^'^^ , 
all with first excited state separations around 150 cm~^. 
However, the Raman spectrum of Ho2Ti207 docs not 
show any clear inelastic light scattering from CF levels 
to compare these calculations to. Although there is some 
weak intensity around ~ 150 cm~^, which compares fa- 
vorably with all of the calculated level schemes, the inten- 
sity of this scattering is insufficient to definitively ascribe 
it to a CF level. 
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FIG. 6: (Color online) The low wavenumber region of the perpendicular polarization Raman spectrum of Tb2Ti207 at 4 K. 
A fitted lorentzian centered at zero (Rayleigh line) is subtracted from the data for clarity. The crystal field levels (see fig. 
[5]) are indicated by the filled lorentzian lineshapes that were fit to the data. A distinction has been made between CF levels 
that have been calculated and observed before (yellow), CF levels that have been calculated before, but so far have not been 
observed (blue), CF levels that have been observed before, but have not been calculated (green) and additional CF levels that 
are observed for the first time using inelastic light scattering (red). 



IV. CONCLUSIONS 

To summarize, several members of the rare earth ti- 
tanates family R2Ti207 were studied using magnetic 
susceptibiUty measurements and polarized inelastic light 
scattering experiments. Lattice excitations were found 
to vary only slightly between crystals with different 
rare earth ions. Temperature dependent measure- 
ments also revealed completely analogous behavior of the 
phononic excitations, except for two anomalous phonons 
in Tb2Ti207, which seem to be coupled to the low energy 
crystal field excitations in that compound. Such crystal 
field excitations were observed in Tb2Ti207 and possi- 
bly also in Dy2Ti207. Only one non-phononic excitation 
is clearly observed in Dy2Ti207, its energy consistent 
with estimates of the first excited crystal field level. For 
Tb2Ti207, all of the the previously determined crystal 
field energy levels were confirmed. Moreover, the result- 
ing energy level diagram was expanded by three newly 
observed CF levels, only one of which has been calculated 
before. Also, one previously reported level was found to 



be both Raman and dipolc active, contradicting its cur- 
rent presumed symmetry. These findings may reflect the 
existence of two inequivalent Tb sites in the low tempera- 
ture structure or a static symmetry reduction of another 
nature, suggesting the recently found structural correla- 
tions to induce a minute static disorder in Tb2Ti2 07 at 
very low temperatures. The new crystal field informa- 
tion may serve to help elucidate the complex theoretical 
enigma of Tb2Ti207. 
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TABLE II: Experimentally determined and literature values of the frequencies (in cm ^) of observed Raman modes in R2Ti2 07 at room temperature and their 
respective symmetry assignments (between brackets). 





This work 




Saha et al^ 


Zhang et al-^ 


Mori et al^ 


Hess et al^ 


Gupta 


et al^ 


Vandenborre et al^ 


R = Gd 


Tb 


Dy 


Ho 


Gd 


Gd 


Gd 


Gd 


Gd [calcd.) 


Gd {obs.) 


Gd 


-104" 


~102" 


~103" 


~105" (-) 


-(-) 


-(-) 


105 (Fig) 


-(-) 


-(-) 


-(-) 


110 (-) 


-128" 


- 


~124" 


-122" (-) 


-(-) 


-(-) 


125 C) 


-(-) 


-(-) 


-(-) 


215 (-) 


205 


209 


212 


214 (F2,) 


215 (F29) 


210 iF2g) 


222 (F29) 


219 {F2g) 


230.6 {F2g) 


225 {F2g) 


225 {F2g) 


260 


256 


269 


297 (F29) 


-(-) 


-(-) 


-(-) 


-(-) 


-(-) 


-(-) 


-(-) 


310 


303 


308 


311 (Fig) 


312 (F29) 


310'' (Eg) 


311 (Eg) 


312 {Eg) 


318.0 {Eg) 


-(-) 


317 {F2g) 


325 


313 


328 


329 (Eg) 


330 {Eg) 


310" {F2g) 


-(-) 


-(-) 


328.2 {F2g) 


317 {F2g) 


347" {Eg) 


450*^ 


450" 


45r 


452" (-) 


-(-) 


470 (Fig) 


452 (-) 


455 [F2g) 


-(-) 


-(-) 


-(-) 


- 


- 


- 


-(-) 


-(-) 


-(-) 


-(-) 


-(-) 


522.2 {F2g) 


-(-) 


-(-) 


517 


518 


519 


522 (Aig) 


518 {A^g/F2g) 


520 (A^g) 


518 (F29) 


519 {A,g) 


526.8 {Atg) 


515 {Alg) 


515 {Aig+F2g) 


554 


557 


550 


562 {F2g) 


547 {F2g) 


570 {F2g) 


585 (F29) 


549 {F2g) 


594.0 {F2g) 


580 (F29) 


580 {F2g) 


~677^ 


~689^ 


~693-^ 


-701-^ (-) 


684 (-) 


-(-) 


-(-) 


680 (F29) 


-(-) 


-(-) 


-(-) 


~703^ 


~706^ 


~720^ 


~724-^ (-) 


701 (-) 


-(-) 


-(-) 


-(-) 


-(-) 


-(-) 


705 (-) 



"These modes are very weak and barely resolved, therefore their exact position is estimated. 

'Mori et al. ascribe this band to trace amounts of Gd203. 

"This is a very weak mode, barely resolved above the noise. Mori et al. ascribe this mode to trace amounts of Ti02. 

''Zhang et al. indicate the lower and higher wavenumber components of this band as the Eg and F2g modes, respectively. 

■^This band has been calculated, rather than observed by Vandenborre et al. 

^Those overlapping modes comprise a weak band, in which the two modes cannot be separately resolved. 



